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The standard theory of electromagnetic cascades onto a photon background predicts a quasi- 
universal shape for the resulting non-thermal photon spectrum. This has been applied to very 
disparate fields, including non-thermal big bang nucleosynthesis (BBN). However, once the energy 
of the injected photons falls below the pair-production threshold the spectral shape is much harder, a 
fact that has been overlooked in past literature. This loophole may have important phenomenological 
consequences, since it generically alters the BBN bounds on non-thermal relics: for instance it allows 
to re-open the possibility of purely electromagnetic solutions to the so-called “cosmological lithium 
problem”, which were thought to be excluded by other cosmological constraints. We show this with 
a proof-of-principle example and a simple particle physics model, compared with previous literature. 

PACS numbers: 95.30.Cq 26.35.-|-c, 14.80.-j 


I. INTRODUCTION 

Electromagnetic cascades, namely the evolution of 
7, particle numbers and energy distribution following 
the injection of a energetic 7 or e in a medium filled with 
radiation, magnetic fields and matter, is one of the physi¬ 
cal processes most frequently encountered in astroparticle 
physics, in domains as disparate as high-energy gamma- 
ray astrophysics, ultra-high-energy cosmic ray propaga¬ 
tion, or the physics of the early universe. In particular, 
the elementary theory of such a cascade onto a photon 
background has been well known since decades, and can 
be shown via a textbook derivation (see Chapter VIII 
in PQ, for instance) to lead to a universal “meta-stable” 
spectrum—attained on timescales much shorter than the 
thermodynamical equilibration scale—of the form: 

! K 0 {— ) for EL < ex , 

2 7 a, 

K 0 for e x < E 1 < e c , 

0 for E > e c . 

In the above expression, K 0 = E 0 e^[ 2 + ln(e c / ex)\~ x is 
a normalization constant enforcing the condition that the 
total energy is equal to the injected electromagnetic en¬ 
ergy, E 0 - the characteristic energy e c = nij/f™ denotes 
the effective threshold for pair-production (e“ ax being 
the highest energy of the photon background onto which 
pairs can be effectively created); ex e c /3 is the maxi¬ 
mum energy of up-scattered inverse Compton (IC) pho¬ 
tons. Natural units with c = ks = 1 are used through¬ 
out. 

A notable application of this formalism concerns the 
possibility of a non-thermal nucleosynthesis phase in the 
early universe (for recent review on this and other as¬ 
pects of primordial nucleosynthesis, or BBN, see mm)- 
The determination of the baryon energy density of the 
universe fib inferred from the CMB acoustic peaks mea¬ 
surements can be used in fact to turn the standard BBN 
into a parameter-free theory. The resulting predictions 


for the deuterium abundance (or 2 H, the most sensitive 
nuclide to D;,) are in remarkable agreement with ob¬ 
servations, providing a tight consistency check for the 
standard cosmological scenario. The 4 He and 3 He yields 
too are, broadly speaking, consistent with this value, al¬ 
though affected by larger uncertainties. The 7 Li predic¬ 
tion, however, is a factor ~ 3 above its determination in 
the atmosphere of metal-poor halo stars. If this is inter¬ 
preted as reflecting a cosmological value—as opposed to 
a post-primordial astrophysical reprocessing, a question 
which is far from settled mm — it requires a non-standard 
BBN mechanism, for which a number of possibilities have 
been explored Ell¬ 
in particular, cosmological solutions based on electro¬ 
magnetic cascades have been proposed in the last decade, 
see for instance [6j. However, typically they do not ap¬ 
pear to be viable [3], as confirmed also in recent inves¬ 
tigations (see for instance Fig. 4 in 7], dealing with 
massive “paraphotons”) due to the fact that whenever 
the cascade is efficient in destroying enough 7 Li, the de¬ 
struction of 2 H is too extreme, and spoils the agreement 
with the CMB observations mentioned above. Actually, 
this tension also affects some non-e.m. non-thermal BBN 
models, see for instance [8]. 

This difficulty can be evaded if one exploits the prop¬ 
erty that ' Be (from which most of ' Li come from for 
the currently preferred value of fib, via late electron cap¬ 
ture decays) has the lowest photodissociation threshold 
among light nuclei, of about 1.59 MeV vs. 2.22 MeV for 
next to most fragile, 2 H. Hence, to avoid any constraint 
from 2 H while being still able to photo-disintegrate some 
'Be, it is sufficient to inject photons with energy 1.6 < 
E 1 /MeV < 2.2, with a “fine-tuned” solution (see e.g. the 
remark in [5] or the discussion in cm). Nonetheless, it 
turns out to be hard or impossible to produce a sizable 
reduction of the final 7 Li yield, while respecting other 
cosmological bounds, such as those coming from extra 
relativistic degrees of freedom (N e g) or spectral distor¬ 
tions of the CMB. A recent concrete example of these 
difficulties has been illustrated in [Sj, which tried such a 
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fine-tuned solution by studying the effects of 0(10) MeV 
sterile neutrino decays. 

In this article, we point out that, depending on the 
epoch, at sufficiently low energies of injection the cas¬ 
cade develops differently and the final spectrum is signif¬ 
icantly altered with respect to Eq. (Qp, which has been 
incorrectly used till recently, see e.g. fO HO] . As a con¬ 
crete application, we show how this re-opens a window 
to a cosmological solution to the 7 Li problem via e.m. 
decays. Additionally, one expects peculiar signatures as¬ 
sociated to such scenarios, which can be probed with cos¬ 
mological observations. We will discuss this both in a 
proof-of-principle example and in the context of a par¬ 
ticle physics model, involving one sterile neutrino. This 
was chosen for its simplicity and to allow for a direct 
comparison with the results of [9], which studied a sim¬ 
ilar model. Further considerations on some additional 
implications of our insight are finally outlined. 


II. E.M. CASCADES AND UNIVERSAL 
NON-THERMAL SPECTRUM 


Our argument is the following: Let us assume that one 
injects photons at some time (or corresponding plasma 
temperature T ) whose energy Eq is below the pair pro¬ 
duction threshold at that epoch, which can be esti¬ 
mated for the CMB plasma to be e c ~ ml/(22 T) ~ 
lOTJ^y MeV [TT]. Note that as long asT < few keV, this 
is compatible with the typical nuclear photo-disintegration 
energies relevant for BBN. It is clear that the spectrum of 
Eq. § cannot stay valid in this regime: there is no pair- 
production cutoff, of course, but even the lower-energy 
part cannot be correctly captured by Eq. ([I]). Unless one 
considers other physical processes for the photon interac¬ 
tions, not included in the derivation of Eq. ([l]), there are 
no non-thermal electrons which can up-scatter CMB pho¬ 
tons! Since the photon interaction probability is much 
smaller below pair production threshold, at leading level 
the injected spectrum below e c stays the same—apart 
for red-shifting, which happens on very long timescales 
with respect to particle photon interactions and hence we 
neglect. Accounting for the finite probability for the pho¬ 
tons to scatter- via 77 , via Compton scattering off the 
background electrons, or via Bethe-Heitler e ± produc¬ 
tion onto background protons and Helium nuclei—one 
does end up with a suppression of the injected spectrum, 
plus a lower energy tail due to downgraded energy 7 ’s as 
well as 7 ’s produced via IC by the secondary e’s. The re¬ 
sulting secondary or tertiary photons, on the other hand, 
are typically at too low-energies to contribute to photo¬ 
dissociations and will be neglected. Within this approxi¬ 
mation, the Boltzmann equation describing the evolution 
of the distribution function / 7 reads: 


0/7^7) 


-r y (E y ,T(t))f^,T(t))+S(E^,t ), ( 2 ) 


where S(E 7 , t) is the source injection term, T is the total 
interaction rate, and we neglected the Hubble expansion 
rate , since interaction rates are much faster and rapidly 
drive / 7 to a quasi-static equilibrium, c> ^'g/~ ^ = 0. Thus, 
we simply have : 


/ 7 s (£ 7 ,f) 


S(E 11 1) 
r7(^7) ’ 


( 3 ) 


where the term S for an exponentially decaying species 
with lifetime Tx and density nx(t), whose total e.m. en¬ 
ergy injected per particle is E 0: can be written as 


5(U 7 ,t) 


n°Cx(l +z(t)) 1 * 3 e t/Tx 

E 0 t x 


Vi(E 1 ) , 


( 4 ) 


with z(t) being the redshift at time t, and the energy 
parameter Qx (conventionally used in the literature) is 
simply defined in terms of the initial comoving density 
of the X particle n x and the actual one of the CMB, 
n°, via n x = n®(/x/E o- A monochromatic emission line 
would then correspond to p 1 (E 1 ) = 6(Ej — Eq). For a 
two body decay X —>■ 7 U into a monochromatic line plus 
another not better specified (quasi)massless particle U, 
one would have E 0 = nix /2, where rrix is the mass of the 
particle. Here, we will be interested in multi-MeV values 
for the mass mx and at temperatures of order few keV 
or lower, hence the thermal broadening is negligible and 
a Dirac delta spectrum as the one above is appropriate. 

The interaction rate T 7 is computed by accounting for: 
i) Compton scattering over thermal electrons, 7 + e t .h —+ 
7 + e, taken from HU; ii) scattering off CMB photons: 
7 + 7 th -> 7 + 7 , for which we follow [12]; iii) Bethe- 
Heitler pair creation : 7 + TV —»• X + e ± , for which we 
use the formulae of [H]. Note that we neglect the small 
effect due to the finite probability for the secondary or 
tertiary photons to induce some dissociations, i.e. once a 
photon interacts it is “lost”. The results that we obtain 
are in this respect slightly conservative, by an amount 
which we estimated to be of the order of a few %. 


III. NON-THERMAL NUCLEOSYNTHESIS 


At temperatures of few keV or lower, the standard 
BBN is over, and the additional nucleosynthesis can be 
simply dealt with as a post-processing of the abundances 
computed in the standard scenario. The non-thermal 
nucleosynthesis due to electromagnetic cascades can be 
described by a system of coupled differential equations of 


1 In the case of the universal spectrum, this hypothesis of quasi¬ 

static equilibrium has been checked (see EH) and is in very good 
agreement with generated spectrum calculated numerically when 

this assumption is not made (see e.g. m or EH). 


dt 
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the type 

( ]Y > _ r°° 

—— =^2Yt dE 1 f 1 {E 1 )a 1+T ^A{E 1 ) 

rj-i J 0 

/»oo 

- hiV / dE 1 f 1 (E 1 )cj 1+ A^p(E 1 ) (5) 

p J 0 

where: Ya = riA_/rib is the ratio of the number den¬ 
sity of the nucleus A to the total baryon number den¬ 
sity rib (this factors out the trivial evolution due to the 
expansion of the universe); (t 1+ t^,a is the photodisso¬ 
ciation cross sections onto the nuclei T into the nucleus 
A , i.e. the production channel for A; a 1 j r A^p is the 
analogous destruction channel (both cross sections are 
actually vanishing below the corresponding thresholds). 
In general one also needs to follow secondary reactions of 
the nuclear byproducts of the photodissociation, which 
can spallate on or fuse with background thermalized tar¬ 
get nuclei (see for instance | 6 ]) but none of that is rele¬ 
vant for the problem at hand. If the injected energy is 
1.59 < Eg/MeV < 2.22, the only open non-thermal BBN 
channel is 7 + 'Be —> 3 He + 4 * He, whose cross-section 2 we 
denote with ay, there are no relevant source terms and 
only one evolving species (since Yy <C 13,4), thus yielding 
for the final (at Zf) to initial (at Zi) abundance ratio 


hr \) = r . 

\Y 7Be (z f )J J E 0 H?t x T(E 0 ,z ) 

( 6 ) 

To obtain Eq. (| 6 |), we transformed Eq. (|5|) into redshift 
space, defining H(z) = 1 + z) 2 as appropriate for a 

Universe dominated by radiation, with H ° = Hoy/Q%, 
Hq and f2/ being the present Hubble expansion rate and 
fractional radiation energy density, respectively. By con¬ 
struction, equating the suppression factor given by the 
RHS of the Eq. ([b]) to ~ 1/3 provides a solution to the 
7 Li problem which is in agreement with all other con¬ 
straints from BBN. In Fig. [TJ the lower band shows for 
each tx the range of ( x corresponding to a depletion 
from 40% to 70%, for the case E 0 = 2MeV. Similar re¬ 
sults would follow by varying Eg by 10 % about this value, 
i.e. provided one is not too close to the reaction thresh¬ 
old. The upper band represents the analogous region if 
we had distributed the same injected energy according 
to the spectrum of Eq. |l]), up to min[e c ,E 0 \. It is clear 
that in the correct treatment a large portion of this region 
survives other cosmological constraints, described below, 
while none survives in the incorrect treatment. 


2 It is worth reporting that the cross-section for this process re¬ 

ported in the appendix of [6] is erroneous. This has already been 

pointed out in [5], which we agree with. The correct formula is 

used in the following. 



FIG. 1: The lower band is the range of abundance parameter 
Cx—, vs. lifetime tx, for which the primordial lithium is depleted 
to 40% to 70% of its standard value, for a monochromatic photon 
injection with energy Eq = 2MeV. The upper band represents the 
analogous region if we had distributed the same injected energy, 
up to Eq = 2 MeV, according to the erroneous spectrum of Eq. (fTji. 
Above the solid blue curve, a change in entropy (and H/,) between 
BBN and CMB time larger than the 2rr error inferred from CMB 
would be obtained. The region to the right of the dot-dashed green 
curve is excluded by current constraints from //-distortions in the 
CMB spectrum El according to the computation of m, while the 
dashed cyan curve illustrates the weaker bounds that would follow 
from the less accurate parameterization of n. The dotted red 
curve is the forecasted sensitivity of the future experiment PIXIE, 
corresponding to |/x| ~ 5 X 10' - 8 Eq]- 


IV. CMB CONSTRAINTS 


We mentioned the baryon abundances Of, inferred from 
CMB and BBN (notably 2 H) probes are consistent within 
errors. This implies that no major injection of entropy 
took place between the BBN time and the CMB epoch, 
otherwise the baryon-to-photon ratio (proportional to 
fib) would have changed, see for instance [T5]. In a 
radiation-dominated Universe, the change in entropy as¬ 
sociated to a release of energy into all e.m. particles 
characterized by parameter Cy-s -e.m. and a lifetime t x 
can be estimated as 


— ~ l n = 2 14 x 10 " 4 — 

S ~ Si ' 10 - 9 GeV 



( 7 ) 


For illustration, in Fig. jT| the solid blue line represent 
the level of entropy release associated to a variation of 
2 a around the best-fit measured value of 1% by Planck, 
A S/S ~ 0.022 pi/ Since the level of injected energy 
needed to solve the lithium problem via a monochromatic 
line is up to two orders of magnitude below the bounds, 
it is clear that this constraint is very weak, but for very 
short lifetimes of the order of 10 4 s. 

Another constraint comes from the level of spectral 
distortions in the CMB. For the relatively short life¬ 
times relevant for the problem, the Compton scatter¬ 
ing is fast enough that energy-redistribution is effective, 
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no y-type distortion survives. On the other hand, pro¬ 
cesses that change the number of photons are relatively 
rare, and a residual distortion of the /x- type is possi¬ 
ble. This has been constrained by COBE-FIRAS to be 
|/x| < 9 x 10" 5 [T7]. The level of spectral distortion pro¬ 
duced by the decay process here has been estimated in the 
past (see for instance m ), but a recent re-evaluation | 18 j 
found significant improvements at short lifetimes, essen¬ 
tially due to a better treatment of the time-dependence 
of the visibility function. The theoretical expectation for 
/x can be written as 

< 8 > 

where the function J is taken from m■ The bound ex¬ 
cludes the region to the right of the dot-dashed, green 
curve in Fig. [I] For comparison, the dashed cyan curve 
reports the much weaker bound that would follow from 
the approximations in EH3- We also checked that the ex¬ 
tra constraint due to extra “dark radiation” parameter¬ 
ized by iV e ff is irrelevant as long as the branching ratio 
in extra relativistic species is not larger than a couple 
of orders of magnitudes with respect to the photon one. 
We thus conclude that there is a significant interval of 
lifetimes (10 4 < tx/s < 10 6 ) and corresponding energy 
injection parameter 10 -3 > (x-^/MeV > 1.3 x 10 -6 for 
which a perfectly viable solution is possible. We remind 
once again that this possibility appeared to be closed due 
to the use of Eq. 0 beyond its regime of applicability. 

One may wonder how realistic such a situation is in a 
concrete particle physics model. Although we refrain here 
from detailed model-building considerations, it is worth 
showing as a proof-of-principle that models realizing the 
mechanism described here while fulfilling the other cos¬ 
mological constraints (as well as laboratory ones) can 
be actually constructed. Let us take the simplest case 
of a sterile Majorana neutrino with mass in the range 
3.2 < M s / MeV < 4.4, mixing with flavour a neutrinos 
via an angle 8 a . We also define 0 2 = The three 

main decay channels of this neutrino are (see e.g. [2T1 and 
refs, therein): 

G 2 M 5 © 2 

• V s ->• 3ix, with rate r„ a _> 3 „ ~ f 92 ^ 3 ; 

• v s —► v a e + e~, with a rate depending on single 8 a ' s; 

• u s —► txy, with a rate i — 256 tt 4 ^ ■ 

The resulting branching ratios for the masses of interest 
and 8 e <C 0 are of the level of 0.9 : 0.1 : 0.01, respectively. 
It is physically more instructive to normalize the abun¬ 
dance of the zx s , ro 3 , in terms of one thermalized neutrino 
(plus antineutrino) flavour species, n 3 . In Fig. [ 2 ] we show 
the corresponding range of parameters in the 0 — n 3 /n 3 
plane, for M s = 4.4MeV, for which the ‘Li problem 
is solved, fulfills cosmological constraints and, provided 
that 9 e <C 0, also laboratory ones [22]. It is worth not¬ 
ing that: i) the entropy release bound is now close to the 



FIG. 2: Constraints for the sterile neutrino model discussed 
in the text. The legend is the same as for case 1). 

region of interest, since the decay mode v s —> n a e + e~, 
which is useless as far as the 'Be dissociation is con¬ 
cerned, dominates the e.m. energy injection, ii) A non- 
negligible fraction of relativistic “dark radiation” is now 
injected, mostly via the dominant decay mode v s — > 3iq 
hence we added the current la sensitivity of Planck to 
N e tf [IS], with A N e g computed similarly to what done 
in [2]- The needed abundance could be obtained in sce¬ 
narios with low reheating temperature j22|- 

V. CONCLUSIONS 

We have discussed the breaking of the universality of 
the photon spectrum in electromagnetic cascades, when 
the energy of the injected photons falls below the pair- 
production threshold. This may be of interest for a 
number of astroparticle applications, but in the spe¬ 
cific case of the cosmological context, this happens when 
E 1 < m 2 /(22T) ~ 10 MeV. We noted that the en¬ 
ergies concerned are of the same level of the binding ener¬ 
gies of light nuclei. This implies a potential large impact 
on non-thermal nucleosynthesis models, notably of elec¬ 
tromagnetic type, but could be also relevant for models 
with late time hadronic cascades. We provided an analyt¬ 
ical estimate of the resulting (much harder) meta-stable 
spectrum of non-thermal photons, and showed that the 
impact is so large that it can potentially re-open the 
possibility of electromagnetic cacade solutions to the so- 
called “lithium problem”, which were thought to be ex¬ 
cluded by other cosmological constraints. We substan¬ 
tiated this point with a proof-of-principle example of a 
photon line injection at ~ 2 MeV from a particle decay, 
satisfying by construction all other BBN constraints but, 
not trivially, also all other cosmological bounds plaguing 
previous attempts. Although we did not indulge into par¬ 
ticle model building, we proved that the right conditions 
can be actually satisfied in a simple scenario involving a 
~ 4 MeV sterile neutrino mostly mixed with v T and/or 
with effective mixing angle 0 ~ 10 -2 . 
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The possibility to find new mechanisms to deplete the 
standard BBN prediction of lithium abundance in a con¬ 
sistent way is probably the most spectacular consequence 
of our investigation. In turn, this could stimulate more 
specific model-building activities. For instance, decays 
of relatively light new neutral fermionic particles A' for 
which the v + 7 channel is the only two body standard 
model channel opened—as it is the case for the light grav- 
itinos in supergravity models—constitute a natural class 
of candidates. Alternatively, one may think of decaying 
scenarios involving a pair of quasi degenerate mass states 
X and Y. which are potentially much heavier than the 
MeV scale. Some of these scenarios may be motivated 
by other astroparticle or particle physics reasons and 
certainly deserve further investigation. We also showed 
how improvements in the determination of /z—type spec¬ 
tral distortions bounds of the CMB might be crucial to 
test these scenarios: testing frameworks for the parti¬ 
cle physics solutions to the lithium problem may thus 
provide additional scientific motivations for future in¬ 
struments like PIXIE :2Q]. Computations of distortions 


corresponding to specific injection histories may also be 
refined: for instance, for short lifetimes relativistic cor¬ 
rections to the double Compton and Compton scatter¬ 
ing may be important to improve the theoretical accu¬ 
racy [531 . 

Finally, from a phenomenological perspective, an ob¬ 
vious spin-off of our work would be to re-compute the 
BBN bounds to electromagnetic decaying particles in 
cases where the universality of the spectrum of Eq. ([I]) 
breaks down. Preliminary results indicate that bounds 
can be easily modified by one order of magnitude. These 
results will be reported in a forthcoming publication. 
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